International Journal of Theoretical Physics manuscript No. 

(will be inserted by the editor) 



Asymptotic behavior of Heun function and its integral 
formalism 

Yoon Seok Choun 



Received: date / Accepted: date 



Abstract Heun function generalizes all well-known special functions such as: Spheroidal 
Wave, Lame, Mathieu, and hypergeometric 2F1, 1^1 and ofi functions. Heun func- 
tions are applicable to diverse areas such as theory of black holes, lattice systems in 
statistical mechanics, solution of the Schrodinger equation of quantum mechanics, 
addition of three quantum spins. 

In this paper, applying three term recurrence formula [1,1 consider asymptotic 
behaviors of Heun function and its integral formalism including all higher terms of 
A„'s. I will show how the power series expansion of Heun functions can be converted 
to closed-form integrals for all cases of infinite series and polynomial. One interesting 
observation resulting from the calculations is the fact that a 2F1 function recurs in 
each of sub-integral forms: the first sub-integral form contains zero term of A' n s, the 
second one contains one term of A„'s, the third one contains two terms of A„'s, etc. 

In section 5, 1 apply my integral formalism of Heun function to "The 192 solu- 
tions of the Heun equation" 11251 . Due to space restriction final equations for all 192 
Heun functions is not included in the paper, but feel free to contact me for the final 
solutions. Section six contain two additional examples using integral forms of Huen 
function. 

This paper is 4th out of 10 in series "Special functions and three term recurrence 
formula (3TRF)". See section 8 for all the papers in the series. Previous paper in 
series deals with the power series expansion in closed forms of Heun function. The 
next paper in the series describes analytically the power series expansion of Mathieu 
function and its integral formalism. 
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1 Introduction 

The Heun function, having three term recurrence relations, are the most outstand- 
ing special functions in among every analytic functions. Due to its complexity Heun 
function was neglected for almost 100 years[4|. According to Whittaker's hypothesis, 
'The Heun function can not be described in form of contour integrals of elementary 
functions even if it is the simplest class of special functions' 

Recently Heun function started to appear in theoretical modern physics. For ex- 
ample the Heun functions come out in the hydrogen-molecule ion 1 17 1, in the Schrodinger 
equation with doubly anharmonic potential[24| (it's solution is the confluent forms 
of Heun function), in the Stark effect lfTBI . in perturbations of the Kerr metric lfT8l[T9l 
l20ll2Tll22l . in crystalline materials! 16], in Collogero-Moser-Sutherland systems[23|, 
etc., just to mention a few. H3III1 Traditionally, we have constructed all physical 
phenomenons by only using two term recursion relation in power series expansion 
until 19th century. However, since modern physics (quantum gravity, SUSY, general 
relativity, etc) come out of the world, we have at least three or four recurrence re- 
lations in power series expansions. Furthermore these type of problems can not be 
reduced to two term recurrence relations by changing independent variables and co- 
efficients. 

In previous paper I showed the exact analytic solution of Heun function for all 
higher terms of A„'s by applying three term recurrence formula[l |; the power series 
expansion of infinite and polynomial cases[ 3]. Now, I consider integral forms of Heun 
function and its asymptotic behavior of it and the boundary condition for the inde- 
pendent variable x. Expressing Heun function in integral forms resulting in a precise 
and simplified transformation of Heun function to other well-known special func- 
tions such as: hypergeometric function, Mathieu function, Lame function, confluent 
forms of Heun function and etc. Also, the orthogonal relations of Heun function can 
be obtained from the integral forms. 

In Ref.[4], Heun's equation is a second-order linear ordinary differential equation 
of the form 



With the condition e = a + — 7— 8 + 1. The parameters play different roles: a 7^ 
is the singularity parameter, a, /3, 7, 8, £ are exponent parameters, q is the accessory 
parameter. Also, a and /3 are identical to each other. The total number of free param- 
eters is six. It has four regular singular points which are 0, 1, a and °° with exponents 
{0, 1 — 7}, {0, 1 — 8}, {0, 1 — e} and {a, j3}. Assume that y(x) has a series expansion 
of the form 




(1) 



y(*) = I 



.n+X 



(2) 



;i=0 
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Plug © into (Q} . 

c n+ \ = A„c„+B n c n -\ ;«>1 (3) 

where, 

{n+X){n-\ + y+e+X+a{n-\+y+X + 8))+q 



A " ' a{n + \+l){n + y+X) 

(n + X)(n + a + j5-8 + X+a(n + 8 + y-l+X))+q 



B n — 



fl(« + l+A)(« + y+A) 

(»-l+A)Q + y+g + £-2 + A) + a/3 _ (»- 1 +A + q)(n- 1 +A + j3) 

a(n + l+A)(n + y+A) ~~ o(n + l+A)(« + y+A) 

(4b) 

ci=A c (4c) 
We have two indicial roots which are Ai =0 and A2 = 1 — y 

2 Asymptotic behavior of the function y(x) and the boundary condition for x 

2. 1 Infinite series 

Now let's test for convergence of infinite series of the analytic function y(x). As n 
goes to infinity, (f4ab and ( f4bT > are 

limA„ = A = -i '- limB„ =B = - (5) 

n»l fl n»l a 

Substitute © into ©. For « = 0, 1 , 2, • • • , it give 



Co 










Ci 


= AC 








c 2 


= (A 2 ^ 


-B)C 






c 3 


= (A 3 ^ 


-2AB)C 




C A 


= (A 4 ^ 


- 3A 2 B + 


B 2 )C 




c 5 


= {A 5 - 


-4A 3 5 + 


3AB 2 )C 


c 6 


= (A 6 - 


-5A 4 B + 


6A 2 B 2 + 


B 3 )C 


c 7 


= (A 7 - 


-6A 5 B + 


10A 3 5 2 - 


f4AB 3 )C 


c 8 


= (A 8 ^ 


-7A 6 B + 


15A 4 fi 2 - 


f 10A 2 5 3 +fl 4 )C 



(6) 



The sequence c„ consists of combinations A and B in ©. First observe the term inside 
parentheses of sequence c„ which does not include any A„'s in ©: c„ with even index 

(cn,C2,C 4 ,- • •)■ 
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(1) Zero term of A's 

CO 

c 2 = Bc 
C4 = B 2 cq 

C 6 = B 3 C (ns 

eg = B 4 co 
c H) = B 5 c 



When a function y(x), analytic at x=0, is expanded in a power series, we write 



y(x) = E (8) 

m=0 



where 



Put© in © putting m = 0. 



n=0 



yo(*)=c £(2fct 2 ) (10) 

n=0 

Observe the terms inside parentheses of sequence c„ which include one term of A„'s 
in (O: c„ with odd index (c\, C3, cs,- • •)• 
(2) One term of A's 

c\ =Ac 
c 3 = 2ABc Q 
c 5 = 3AB 2 c 
Cl =4AB 3 c 
c 9 = 5AB 4 c 



Put d2~TT > in © putting m = 1. 

^W^coAxf^l^x 2 )" (12) 
«=o 1 • 

Observe the terms inside parentheses of sequence c„ which include two terms of A„'s 
in (|6]i: c„ with even index (C2, C4, ce,- ■ ■ )■ 
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(3) Two terms of A's 

c 2 = A 2 c Q 
c 4 = 3A 2 Bc Q 
c 6 = 6A 2 B 2 c 

c 8 = \QA 2 B i c ( 13 ) 
cio = l5A 2 B 4 c 



Put (O in (O putting m = 2. 

y 2 ( X )^c (A X ) 2 i ^ +1 ^ + 2 \ B X 2 Y (14) 
«=o L - 

By repeating this process for all higher terms of A's, we can obtain every y m (x) terms 
where m > 2. Substitute ( TTOb . | fl2l . ( fT4b and including all y„,(:ic) terms where m > 2 
into ©. 

y(x) = £ c„x n+;i =yoW +yi(x) + yi{x) +y 3 (x) +■■■ 
,1=0 

= Y Y ( n + m ) ] j?f> where c = l,x = Bx 2 and y = Ax (15) 

(fTBI l is simply 

= 1 ; ^ where + < 1 (16) 

1 - (x+y) 

( fT6l l is geometric series. Put x = Bx 2 and y — Ax into the condition of convergence of 

CE6}. 

(1+fl) 2 / l+fl\ 2 (1+fl) 2 

' -a<\x — <- — -r^+a (17) 



4 V 2 J 4 

According to 1171 . I have a boundary condition of x for the infinite series of Heun 
function(Heun function) in which is following the way. 

(A) As a=0 

no solution (18a) 

(B) As a=l 

1-V2<*<1 and 1<*<1 + V5 (18b) 

(C) As < a < 1 



(l+o)- Va 2 + 6a + l , , (l+a) + Va 2 +6a + l 

<x<a and 1 < x < - (18c) 

2 2 

(D) As a > 1 
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[l+a)-^a 2 + 6a + l , , (1 +a) + \/a 2 + 6a + 1 /10J ^ 
<x<l and a<x<- (18d) 



2 

(E) Asa = -3-2V5 



-3-2V2<x<-l-V2 and - 1 - V2 < x < 1 (18e) 

(F) Asa = -3 + 2V5 

-3+2V2<x< -1 + V2 and - 1 + V2 < x < 1 (18f) 

(G) As -3-2V2 < a < — 3 +2-\/2 

a <x < 1 (18g) 

(H) As -3 + 2\/2 < a < and a < -3 - 2\/2 

(1+a)- Va 2 +6a + l , (l+a) + Vfl 2 +6fl + l , /10lx 

a<x<- and - <x<l (18h) 

2 2 

2.2 The Case of a w — 1 

Let assume that a is approximately close to — 1 . But a ^ — 1 in (0. Then A„ terms are 
negligible. (01 is approximately 

c„ +1 ssB„c„_ 1 (19a) 

And, 

\imB„fzB = - (19b) 

n>l a 

We can classify c n as to even and odd terms from plugging (119b) into (119a) . 



CO 




c\ 




C2 


= Bc 


C3 


= Bci 


C4 


= B 2 c 


C5 


= B 2 c { 


C6 


= B 3 c 


Cl 


= B 3 Cl 


C8 


= B 4 c 


c 9 


= B A Cy 



(20) 



When a function y(x), analytic at x=0, is expanded in a power series x=0 by using 
(120) . we write 

limy(x)=c £ (Bx 2 )"+cix£ (gx 2 )" = c - _ +ci - _ fo2 where|fix 2 |<l 

n W =0 77=0 ^ ^ 

(21) 
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Substitute ( 1 19bb into (f2TT> . And for simplicity, let say Co = c\ = 1 into it. 

1 +x 



lim y(x) 
«>i v ' 



1 



±* 2 



where a ^ 



The condition of convergence of x is 



— x 

a 



< 1 



(22) 



(23) 



2.3 The Case ofa^-lora^; — 1 

Let assume that a is much greater than 1 or is much less than — 1 . Then B n terms are 
negligible. (0) is approximately 



And, 



limA,, ss A : 



ci =A c « Ac 



(24a) 
(24b) 



Substitute ( I24bt into ( I24ab . For n = 0, 1 ,2, • • • , it give 

Co 

d =AC 
C 2 =A 2 C 
C 3 =A 3 C 
C 4 =A 4 C 



(25) 



When a function y(x), analytic at x=0, is expanded in a power series x=0 by using 
), I write 



limy(x) =c V (Ax) n 



1 



I -Ax 



where \Ax\ < 1 and cq = 1 



Substitute J24bb into (l26b . 



lim y(x) 
«>i v ' 



l-H> 



d27T i is binomial series. The condition of convergence of x is 

(1+fl) 



< 1 



(26) 



(27) 



(28) 



We can obtain accurate numerical values of Heun functions using machine calcula- 
tion from the above all asymptotic cases. Also, we might be possible to obtain ana- 
lytic power series expansions in closed forms of all 192 local solutions of the Heun 
equation analytically 1 25 1. 
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3 Integral Formalism 

3.1 Polynomial in which makes B n term terminated 

3.1.1 The case of a = —2a; — i — X and =/= —2/3; — i — X where i, a,, j3, = 0, 1,2, • • • 

Now let's investigate the integral formalism for the polynomial case of B n term ter- 
minated at certain eigenvalue. There is a generalized hypergeometric function which 
is: In this paper Pochhammer symbol (x)„ is used to represent the rising factorial: 

{X)„ - rjjp 

i/=ij_l ( C^/);'/-! ( f + 3 + + 2 + f + 2 + 3 + fO'i 



'/-I 



(^t+l+D-H^+l-l+l+l)- 1 ^;! 



By using integral form of beta function, 

B + ^ + + 1) = £ dt, ^-^(l-t.y (30a) 

B(/,-i + ^-^ + | + |,j+l)=/ o 1 ^/"; , - 1 + ^ i+l7+4 (l-"/V (30b) 

Substitute (f30ab and (l30bl into (|29]l. And divide i + | + § ) (//_ 1 + 3- 3 + 2 + 3") 
into //. 

1 



(/ / _ 1 + i + |)(/ / _ 1 + i-± + f + §) 
x ^ (~ a /)»/(f + 3 + ?)'>(! + I + (3 + 2 + 3 + 



A/ tf 1+T / dui uf 2+5 +i (ztiui) 



2 T 2 T 2/'l-n' T 2 1 2^^2 T 2 T 2 n 2 /</ 
1 3 A 



'/-I 



x £ fizlZ^fe_l±i±i±i)z [z(1 _„)(, „„ 

;=0 



The integral form of hypergeometric function is 



l r(i-o)r(y) 



27Ti r(y-a 
whereRe(7-a) >0 



^rfv, (-vij^^l -v/) 5 '-"- 1 ^ -zv/)~!32) 
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replaced a, j3, yandz by z*,_i - a/, z',_i + | + | + 1, 1 andz(l -f/)(l -m/) in <[32j» 



£ft zl ^Mti±i±i±i)i [z(1 _ ( , )(1 _„, ) H 

^•/^I(l-I) ai (l- Z v /( 1-,)(1- M; ))-^ 
vi 1 



(v,-l)l-zv,(l -//)(! -ii,) 



(33) 



Substitute ([33]) into ( TO . 



1 



j4 ( Oij);, ( 2 + 2 + 2 ),■,(! + 2+2 )'/-! (2 + 2 + 2 + 2 )'/-! , 

k=H-i (~ a i)ii-i (2 + 1 + (1 + 1 + 1)11(2 + 2 + 2 + §)'; 



"1 / ,,i /-l ; 3,7, a 1 /• 1 / 1 \ °' 

= dt.tr 1 / t/M/ U.f~ — — <E C/V; — ( 1 — ■ » 



Jo 2ni J v; V v, 



10 



Yoon Seok Choun 



In Ref . [ 3 ] , the general expression of power series of Heun function for polynomial in 
which B n term terminated where a = —2a, — i — X and j3 ^ — 2/3; — i — X is 



(-«o)f (f + |)-o j 

.1 + 2")'o(2 + 2 + 2')'0 

|){'o+2TTT^)(- 2o: o + i 8 - g + a ( g + y- 1 + A ))} + 2Tfey 
0'o + 2 + |)0'o + | + |) 

x (-«0) i0 (§ + |),Q £ f + f + |), L (| + |), Q (1 + | + |) ;o " 

(1 + j)k(l + 2 + t)'o i'i=«d I (~ a i)io(j + f + §)>'o(| + + 3 + 



n=0 






f «0 








I '0=0 1 




O'o + 


+ E 


<0=0 



£ f g (r + |){«o+2(TT5j(-2ab + j3-5+fl(5 + y-l+A))} 



2(1- 



iS ( r^O ('0 + 3 + |)('0 + I + |) 




x (-«b)f (f + |)io 

(1 + 2") 1 'o(2 + 2 + 3") ! 'o 

g (4 + § + |){^+2 (T ^-y(-2^ + j3-g + a (g + y+A + fc-l))} + in ^ y 
*=i l«*=«*-i (4 + 1 + 2 + |))('i + | + 2 + l) 

+ f + +_| + gOjfc— i (; + 1 + 1 + j)'k-i \ 
(— 0£*)i t _i (| + f + f (1 + 1 + §)tt(i + I + 2 + 3f)<* J 
x y (~ a «)» n (| + 2 + l)'^ 1 + 1 + + f + 1 + fk-i ^nl^nl (35) 

'«=i„-i ( _a «)( n -i (§ + f + §)/„-] (1 + f + t);„ (5 + 3 + 3 + t)'/i J J 



X 



where 



(36) 



and 



a = —2a, — i — X as j = 0, 1,2,--- and a; = 0,1,2, 
a, < a 7 - only if i < j where i, j = 0, 1 , 2, • • ■ 



(37) 
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Substitute (l34t into (l35l l where / = 1,2,3, and the integral formalism of Heun 
function for polynomial in which B„ term terminated is 

y( x ) = 

n=0 

-Jiff (-"0k(f + |)^0 J„ 
I i =0 [ l T 2 jio\2 ~T 2 " t " 2 )'0 



in 

«=1 I ir=0 



n-1 / /-l i 



du n h U 1 



n-k U n _ k 



(n-k-3+y+X) 



^—ldv n - k -^—[\ — ] (1 - V„_ jt+ i,„v„_i(l -f„_ fc )(l -«„_*)) ,( " k+ P +X) 

2%i J v„_ k \ v n -kj 



+ 



{-2a n - k -i+p-8 + a(S + r+n-k-2 + X)) 



2(1 +«; 
£ 



2(1 + fl ) 



( =0 (1 + 2)10(2 + 2 + 2 )'o 



(38) 



where 



(39) 



(vi-i) r^riN 

Z only if i > 7 

Put c = 1 as A=0 and c = a~^ 1_r ) as A = 1 - 7 in ([38]). Also, apply ([32]) in it. Then, 
we obtain two independent solutions of Heun equation. The solution is the following 



ways. 

(I) As X = 



1 



; 1 / \ 1 (1 +Cl) l ■> 

y(x) = HF a . p { CCj = - - {a+j) \ j=w „:,T] = — x;z = -x L 



P. 1 , r. 



n-\ / .1 



n=l ^ ifc=0 

a n-k 



-k 



du n _ k u n 



\(n-k-i+y) 



If I / I \ a n-k _ 1 

— idv n ^ k 1 (l- V„_jfc + i^V„_fc(l-^_jfc)(l -U„-k)) - 



1 



2(1 + a 



■(-2a n _ jfe _i + j3-S+a(5 + y+n-fc-2)) 



2(l+a) 



:2Fl (_«,,| ; ^ + | ; V,„)}n» 



(40) 
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y(x) = HS ajtP ccj = --(a + 1 - Y+j)\ j= , 



;'=o,i,2, 

2^1-00,- + ---,---^ 



(1+fl) 1 2 ^ 

; r\ — x;z = -x 

a a , 



n=l I k=0 



(n-k-l-y) 
k 



> /"^l \(n-k-2) 1 / 1 

/„ du n - k u n _ k j-fdv n _ k - 

X (l-^n-k+\,nVn-k{'L-t„- k ){l-U n _ k )y' l{n ~ k+l+ P~ y) 



■(-2a n _£_ 1 +j3-S+a(S+/i-fc-l)) 



2(1 + a 
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2(l+o) 



(l40b is the integral formalism of the first kind of independent solution of Heun func- 
tion for the polynomial as a = —2aj — j where j, 0Cj — 0, 1,2, • • • . And (l4lT i is the 
integral formalism of the second kind of independent solution of Heun function for 
the polynomial as a = — 2a j — j — 1 + 7 where j = 0, 1,2, • • • . 
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3.1.2 The case of (X — —2a,- — i — X and /3 = —2)3; — i — X only if a,- < j3; where 
i,CCi,Pi = 0,1, 

Putj3 =-2fy-i-X where z = 0, 1,2,- •• in®. 



n=0 



A £ (-*M-ft)> 

[ i =0 (! + 2 )'o(2 + 2 + 2 )«0 



H=l I £=0 







n—k 



au n _ k U n _ k 



If 1 / 1 \ a " _ * / 
x— <f> dv n _ k 1 (1 - V„_t + i,„v n _^(l -f„_i)(l -M n _jt)) 

27n j v n _* v v„_ fe y 



-Jt-l+A,) 



+ 



2(1+ a 
2(1 



■(-2a„_ 1 _ A: -2j3„_ 1 _i-5-n + l + Ar-A + a(5 + y+n-A:-2 + A)) 

(42) 



4_1\ 3 (-<*>),•„(-&>)/„ V , 1-1 



Put co= 1 as A=0 and co = a zC 1 as A = 1 — 7 in d42b . Then, we obtain two inde- 
pendent solutions of Heun equation. The solution is the following ways. 
(I) As A = 



y {x) = HF ajPj f ccj = -~(a+j),Pj = hp • ./i , 0) i,2,- ;TJ = -^-^4*;z = -x 2 



2F1 ( -ao,-j3o;i + |;z) +£ J f[ ( / dt n _ k tf" k k 2) [ du^. k u* 



2 2 



n=\ {k=0 



1 /" 1 / 1 \ 

:—&dv n _ k 1 (l-V„ 

2th J v„_ k V v n _fty 



fc-i) 



-fe+i^v„_jfe(l — f„_i)(l —u n _ k )) 



1 



2(1 + a 



■(-2a„_l_t-2j3„_i_i-5-n + l + Ar + a(5 + y+n-fe-2)) 



2 



2(1 +a) 



• :F, ( -ao,-j3 ;i + |;Vi,„) ^>j" 
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(II) As A = 1 - 7 



y(x) = HS ajiPj (aj = ~(a + l - Y+j),Pj = -\w + 1 - r+j)\ J=0AX ... 



(1+fl) 1 2 

;?7 = — -x;z = -xr 

a a 



( — ceo,— j8 ; | — 



B-l / -1 



in 

,i=l I k=0 







£(n-*-2) 
UUn— k u n —k 



1 



2m J v n - k \ v„_£ 



-i(n-t-7) 



\{n-k-y) 



^n-k.n^,, 



2(1 +a) 



{-2a n ^ k -2^ n ^ k -8 + y-n + k + a{8 + n-k-\)) 
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2(1 +a) 



(44) 



(1431 is the integral formalism of the first kind of independent solution of Heun func- 
tion for the polynomial as a = —2dj — j and /3 = — 2j3y — j only if a, < /3y where 
j,OCj,l3j = 0, 1,2, • • • . And (l44i i is the integral formalism of the second kind of inde- 
pendent solution of Heun function for the polynomial as a = —2dj — j—l+Y and 
j3 = -2/9/ - j - 1 + 7 only if a, < where j, o y , j3y = 0, 1 , 2, • • • . 



3.2 Infinite series 



For infinite series, replace the finite summation with an interval [0, Oo] by infinite 
summation with an interval [0,°°] in d38l l. Also replace a,- by — |(a + / + A) on it 
where i = 0, 1 , 2, • • • . Then we obtain the integral formalism of infinite series of func- 
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tion y(x) 



y{x) = Iy»M 

n=0 



I 



f« + Av (£ _i_ Ay 

v 2 ' 2 /'(A 2 ^ 2 JiQ 
10=0 (1 + 2)'o(2 + 2 + 2 )'0 



in 

n=l I Ar=0 



/ /-l 1 



_(n-k-2+X) 
n-k l „-k 



i(n-*:-3+r+A) 



1 -■ 



1 



1 I,/ 1 
2niJ v n _ k \ v n -k 



x (l- V„-jfc + i )B v„_jfc(l-f„_jfc)(l-M„_fc)) 2( " 



(a + j3-S+n- 1 -& + A+a(S + y+n-fc-2 + A)) 



2(1+ a 

(f + l)io(f + f)«0 ^i'o 



9 



2(1 + fl ) 



i =0 (1 + 2 )'0 ( 2 + 2 + 2 )'0 

Put co= 1 as A=0 and cq — cT^ l ~^ as A = 1 — y in ( |451 ). Then, we obtain two inde- 
pendent solutions of Heun equation. The solution is the following ways. 
(I) As A = 



(45) 



y{x) = HF afi \n 

a j5 1 7 



1 



(1+fl) i 2 

x;z = -x 



= 2^1 



2' 2 2 2 
1 



in 

«=1 I k=0 



{(n-k-2) 
dt n -k t n -k 



^(n-k-3+r) 



x <p dv n _ k 

2%iJ v n _ k \ v„_ k 



(1 - V„_ fc+1 ,„V„_ <: (l -/•„-*)(! -Mb-*)) 2( " 



1 



2(1 + a 



(a + ]8 -5 + re- 1 -£ + a(c5 + y+n-fc-2)) 



2(l+o) 



(46) 
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(II) As X = 1 - 7 
y(x) = HS a , p ( rj 



(1+fl) 1 2 

•- -X\Z = —XT 

a a 



2F1 



2 2 2'2 2 2'2 2 



00 f B-l 

in 



x !-■ 



1 I fc=Q 

1 \ -2(»-*+i+«-r) 



n-jfc-2) 1 



2m 



V„-jfc 



1 — ^Ji-jfc+l,rtV n -Jfc(l — ^«-*) (1 —«»-*)) 2 ' + + ^ ^ 



y („-*- y )/ y x y („ , 

n—k,n V " /f -" vi n _fc 5 „ y n—k,n 



-k-y) 



n K,n w) n _ k 



2(1 + fl ) 



(a + J3 - 5 - 7+ n - Ar + a(5 + n - k - 1 )) 



q 



2(1+ a) 



^(f4-i r .M-I=i-I^)M 

d46l l is the integral formalism of the first kind of independent solution of Heun func- 
tion for the infinite series. And (|47T i is the integral formalism of the second kind of 
independent solution of Heun function for the infinite series. 



(47) 



4 Confluent forms of Heun's differential equation 

There are four kinds of confluent forms of Heun equation. [16. 24 28 29 30] We can 
derive these four confluent forms from Heun equation by combining two or more 
regular singularities to each other to take form an irregular singularity. Its process, 
converting Heun equation to other confluent forms, is similar to deriving of confluent 
hypergeometric equation from the hypergeometric equation. First, Confluent Heun 
Equation has regular singularities at x — and 1, and an irregular singularity of rank 
1 at x = 00 as following. 

dZy ■ ( y • 8 '.e]^ + ^- y = (48) 



dx 2 \x x—1 J dx x(x— 1) 

Some examples of the confluent Heun equation are Mathieu functions! 10|, spheroidal 
wave functions lfTTIl . and Coulomb spheroidal functions! 12 1 . 

Doubly-Confluent Heun Equation has irregular singularities at x = and °°, each 
of rank 1 as following. 

d 2 y ( 8 7 \ dy ax-q „ ^ a „ s 

+ I + 1 )-r + — T^y = o ( 4 9) 
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For example, Doubly-Confluent Heun Equation appears in the massive Klein-Gordon 
field on the Kerr spacetime[26]. 

Biconfluent Heun Equation has a regular singularity at x = 0, and an irregular 
singularity at °° of rank 2 as following. 

3-3+ [- + 0+x J— H y = (50) 

This equation is the special case of Grand Confluent Hypergeometric equation as we 
defined before. [2 1 

Triconfluent Heun Equation has one singularity, an irregular singularity of rank 3 
at x = °° as following. 

^+x(Y+x)^- + (ax-q)y = (51) 

As we investigate the Dirac equation in confining potentials ll27l . triconfluent Heun 
equation makes an appearance. 

(I48l>-(f5l1l are confluent forms of Heuns differential equation as we know. We can 
obtain the analytic solutions of these confluent forms of Heun function by replacing 
independent variable x and changing coefficients. Or we are able to have power se- 
ries expansion, integral forms and generation functions of these four second ordinary 
differential equations by using the three term recurrence formula directly. [ 1 1 : if the 
time is permitted, I will publish these four confluent forms of Heun equations. 



5 Integral formalism of 192 Heun functions 

1. A machine-generated list of 192 (isomorphic to the Coxeter group of the Cox- 
eter diagram Da) local solutions of the Heun equation was obtained by Robert S. 
Maier(2007) 11251 . We can obtain integral forms in closed form of all 192 local solu- 
tions of the Heun equation analytically by using three term recurrence formula |Q]; 
the singularity parameter a ^ decides various ranges of independent variable x ac- 
cording to asymptotic behaviors of Heun function. For example, one of the 192 local 
solution of Heun function in Table 2 ll25ll is 

(1 -xy- s Hl(a,q-(8-l)ya;l5-8 + l,a-8 + l,r,2~S;x) (52) 

Replacing coefficients q, a, j3, and 8 by q — (8 — l)ya, — 8 + 1, a — 8 + 1 and 
2-8 into (gOj!, dHJ, (M . (|46T >. (|47I >. we obtain integral forms in closed forms 

of 63. 



6 Additional examples of Heun function in Schrodinger equation and chemistry 

2. We can apply an integral formalism and power series expansion of Heun func- 
tions in many modern physical areas. For example, the Heun functions appear in 
the solution of Schrodinger equation to the quadratic potentials with inverse even 
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powers of two, four and six. |[T3l The solution of the Schrodinger equation to sym- 
metric double Morse potential also need these function. [14] Also, in "The stark effect 
from the point of view of Schrodinger quantum theory [ 15 1, the author considers the 
Schrodinger equation for the hydrogen atom in a constant electric field of magnitude 
E in the z direction. The Schrodinger equation results into two separated equations by 
using parabolic coordinates (see (7), (10) in Ref. lfT5l ). These two equations are of the 
Biconfluent Heun form. Biconfiuent Heun equation can be obtained from Heun equa- 
tion by replacing independent variables and changing coefficients. And as we put the 
new variables and coefficients into integral forms of Heun function on the above for 
the case of polynomials and infinite series, we might be possible to construct power 
series expansions and integral forms in closed forms of Biconfluent Heun function. 
After then, it might be possible to obtain specific eigenvalues for the entire region of r 
by using the power series expansion of Biconfiuent Heun equation. Using the integral 
forms of Biconfiuent Heun equation, it might be possible to construct the normalized 
wave functions and expectation values of any physical quantity as we want. 

3. In "The ionised hydrogen molecule! 17 1, the author consider the hydrogen- 
molecule ion or dihydrogen cation in the Born-Oppenheimer approximation. He 
obtains two individually Confluent Heun equations using the prolate spheroidal coor- 
dinates (see (1), (2) in Ref.| 17]). By replacing independent variables and coefficients 
in Heun equation, we can construct Confluent Heun equation. We might be possi- 
ble to build power series expansions and integral forms in closed forms of Confluent 
Heun function putting the new variables and coefficients into integral forms of Heun 
function on the above for the case of polynomials and infinite series. In general, most 
of wave-functions in physics are quantized with specific eigenvalues. So all solutions 
on the above examples might be quantized with certain eigenvalues. It means that 
its analytic wave-functions have polynomial expansions. And there are infinite num- 
bers of eigenvalues surprisingly because of its three term recurrence form[l |. Also, 
we can transform representations in the form of integrals in Heun function to other 
well-known special functions in an easy way analytically. Because as we see integral 
forms of Heun function, these function include iF\ Hypergeometric function in itself 
on gO]), tED, 63), 61, g6]), g7]). 

7 Summary 

In my previous paper I show the power series expansion in closed forms of Heun 
function (infinite series and polynomial) including all higher terms of A„'s. In this 
paper I derived the integral formalism of Heun function and its asymptotic behaviors 
including all higher terms of A„'s; applying three term recurrence formula flT|. 

As we see the power series expansions of Heun function for all cases of infi- 
nite series and polynomial, denominators and numerators in all B„ terms arise with 
Pochhammer symbol: the meaning of this is that the analytic solutions of Heun func- 
tion can be described as hypergoemetric functions in a strict mathematical way. We 
can express representations in closed form integrals in an easy way since we have 
power series expansions with Pochhammer symbols in numerators and denomina- 
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tors. We can transform Heun function into all other well-known special functions 
with two recursive coefficients because a 2F1 function recurs in each of sub-integral 
forms of Heun function. 

Since we get the integral forms of power series expansions in Heun function, 
we are able to obtain generating functions of it. The generating functions are really 
helpful in order to derive orthogonal relations, recursion relations and expectation 
values of physical quantities. 



8 Series "Special functions and three term recurrence formula (3TRF)" 

This paper is 4th out of 10. 

1 . "Approximative solution of the spin free Hamiltonian involving only scalar po- 
tential for the q — q system" ll3D - In order to solve the spin-free Hamiltonian with 
light quark masses we are led to develop a totally new kind of special function the- 
ory in mathematics that generalize all existing theories of confluent hypergeometric 
types. We call it the Grand Confluent Hypergeometric Function. Our new solution 
produces previously unknown extra hidden quantum numbers relevant for descrip- 
tion of supersymmetry and for generating new mass formulas. 

2. "Generalization of the three-term recurrence formula and its applications" 11321 
- Generalize three term recurrence formula in linear differential equation. Obtain the 
exact solution of the three term recurrence for polynomials and infinite series. 

3. "The analytic solution for the power series expansion of Heun function" Q - 
Apply three term recurrence formula to the power series expansion in closed forms 
of Heun function (infinite series and polynomials) including all higher terms of A„s. 

4. "Asymptotic behavior of Heun function and its integral formalism", 041 - Ap- 
ply three term recurrence formula, derive the integral formalism, and analyze the 
asymptotic behavior of Heun function (including all higher terms of A„s). 

5. "The power series expansion of Mathieu function and its integral formalism", 
ll35l - Apply three term recurrence formula, analyze the power series expansion of 
Mathieu function and its integral forms. 

6. "Lame equation in the algebraic form" 11361 - Applying three term recurrence 
formula, analyze the power series expansion of Lame function in the algebraic form 
and its integral forms. 

7. "Power series and integral forms of Lame equation in the Weierstrass's form 
and its asymptotic behaviors" [37] - Applying three term recurrence formula, derive 
the power series expansion of Lame function in the Weierstrass's form and its integral 
forms. 
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8. "The generating functions of Lame equation in the Weierstrass's form" [38 1 - 
Derive the generating functions of Lame function in the Weierstrass's form (including 
all higher terms of A„'s). Apply integral forms of Lame functions in the Weierstrass's 
form. 

9. "Analytic solution for grand confluent hypergeometric function" ll39l - Apply 
three term recurrence formula, and formulate the exact analytic solution of grand 
confluent hypergeometric function (including all higher terms of A„'s). Replacing ji 
and E(0 by 1 and — q, transforms the grand confluent hypergeometric function into 
Biconfluent Heun function. 

10. "The integral formalism and the generating function of grand confluent hy- 
pergeometric function" Il40l - Apply three term recurrence formula, and construct 
an integral formalism and a generating function of grand confluent hypergeometric 
function (including all higher terms of A„'s). 



Acknowledgment 

I thank Bogdan Nicolescu. The discussions I had with him on number theory was of 
great joy. 



References 

1. Choun, Y.S.: Generalization of the three-term recurrence formula and its applications. 
larXiv:1303 .0806 

2. Choun, Y.S. and Catto, S.: Approximative solution of the spin free Hamiltonian involving only scalar 
potential for the q — q system. arXiv: 1 302.7309 

3. Choun, Y.S. : The analytic solution for the power series expansion of Heun function. arXiv: 1303.0830 

4. Heun, K.: Zur Theorie der Riemann'schen Functionen zweiter Ordnung mit vier Verzweigungspunk- 
ten. Mathematische Annalen 33: 161(1889). 

5. M. Hortacsu: Heun Functions and their uses in Physics. arXiv:l 101.0471 

6. Birkandan, T., Hortacsu, M. : Examples of Heun and Mathieu functions as solutions of wave equations 
in curved spaces. J. Phys. A: Math. Theor. 40, 1105-1116(2007). 

7. Maier, R.S.: The 192 solutions of the Heun equation. Math. Comp. 33, 811-843(2007). 

8. Suzuki, H, Takasugi, E., Umetsu, H: Analytic solution of Teukolsky Equation in Kerr-de Sitter and 
Kerr-Newman-de Sitter Geometries. Prog. Theor. Phys. 102, 253-272(1999). 

9. Suzuki, H, Takasugi, E., Umetsu, H: Perturbations of Kerr-de Sitter Black Hole and Heun's Equation. 
Prog. Theor. Phys. 100, 491-505(1998). 

10. Birkandan, T, Hortacsu, M.: Dirac equation in the background of the Nutku helicoid metric. J. Phys. 
A: Math. Theor. 48, 092301(2007). 

11. Kokkorakis, G. C. and Roumeliotis, J. A.: Electromagnetic eigenfrequencies in a spheroidal cavity 
(calculation by spheroidal eigenvectors). J. Electromagn. Waves Appl. 12(12), 16011624(1998). 

12. Meixner, J., Schfke, F. W. and Wolf, G.: Mathieu Functions and Spheroidal Functions and Their 
Mathematical Foundations: Further Studies, Lecture Notes in Mathematics, Vol. 837, (Springer- 
Verlag, Berlin-New York, 1980). 

13. Figueiredo, B.D.B.: Inces limits for confluent and doubleconfluent Heun equations. J. Math. Phys. 46, 
113503(2005). 

14. Figueiredo, B.D.B.: Generalized spheroidal wave equation and limiting cases. J. Math. Phys. 48, 
013503(2007). 



Asymptotic behavior of Heun function and its integral formalism 



21 



15. Epstein, P.S.: The stark effect from the point of view of Schrodinger quantum theory. Phys. Rev. 2, 
695(1926). 

16. Slavyanov, S.Y. and Lay, W.: Special Functions, A Unified Theory Based on Singularities, (Oxford 
University Press, 2000). 

17. Wilson, A.H.: A generalized spheroidal wave equation, and The ionized hydrogen molecule. Proc. 
Roy. soc. London A 118, 617-647(1928). 

18. Teukolsky, S.A.: Perturbations of a Rotating Black Hole I. Fundamental Equations for Gravitational, 
Electromagnetic, and Neutrino-Field Perturbations. Astroph. J. 185, 635(1973). 

19. Leaver, E.W.: An Analytic Representation for the Quasi-Normal Modes of Kerr Black Holes. Proc. 
Roy. soc. London A 402, 285(1985). 

20. Batic, D., Schmid, H. and Winklmeier, M.: The generalized Heun equation in QFT in curved space- 
times. J. Phys. A 116, 517(2006). 

21. Batic, D. and Schmid, H: Heun equation, Teukolsky equation, and type-D metrics. J. Math. Phys. 48, 
042502(2007). 

22. Batic, D. and Sandoval, M.: The hypergeneralized Heun equation in QFT in curved space-times. 
Central Europ. J. Phys. 8, 490(2010). 

23. Takemura, K.: The Heun equation and the Calogero-Moser-Sutherland system I: the Bethe Ansatz 
method. Commun. Math. Phys. 235, 467-494(2003). 

24. Ronveaux, A.: Hems Differential Equations, (Oxford University Press, 1995). 

25. Maier, R.S.: The 192 solutions of the Heun equation. Math. Comp. 33, 81 1-843(2007). 

26. Dolan, S.R.: Instability of the massive Klein-Gordon field on the Kerr spacetime. Phys. Rev. D 76, 
084001(2007). 

27. Giachetti, R. and Sorace E.: States of the Dirac Equation in Confining Potentials. Phys. Rev. Lett. 101, 
190401(2008). 

28. Lay, W. and Slavyanov, S. Y.: The central two-point connection problem for the Heun class of ODEs. 
J. Phys. A 31(18), 42494261(1998). 

29. Wolf, G.: On the central connection problem for the double confluent Heun equation. Math. Nachr. 
195, 267276(1998). 

30. Bhring, W.: The double confluent Heun equation: Characteristic exponent and connection formulae. 
Methods A ppl. Anal. 1(3) , 348370(1994). 

31. Link to arXiv: 1302.7309 

32. Link to arXiv: 1303.0806 

33. Link to arXiv: 1303.0830 

34. Link to arXiv: 1303.0876 

35. Link to arXiv: 1303.0820 

36. Link to arXiv: 1303.0873 

37. Link to arXiv: 1303.0878 

38. Link to arXiv: 1303.0879 

39. LinktoarXiv:1303.0813 

40. LinktoarXiv:1303.0819 



